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m ■ Abstract 
(N 

^ ■ A complete version of the Moscow A^A^ potential model is presented. The 

■ excellent description for all essential partial waves has been found in the en- 

■ ergy range - 350 MeV. The one-channel version of the model includes the 

■ orthogonality condition to most symmetric six-quark states in all lowest par- 



tial waves and thus, from this point of view, the model generalizes the well 
known Saito's orthogonality condition model (OCM) for the baryon-baryon 
interaction case. The specific features of the presented model which distin- 



^ . guish it from many conventional force models are discussed in details. One 

of them is a specific tensor mixing between nodal and nodeless wavefunctions 
which results in very reasonable values of the OPE cut-off parameter A = 0.78 
GeV and the vrA^A^-coupling constant value = 0.075 in nice agreement with 
modern trends. The model, in case of its confirmation in precise few-nucleon 
calculations, can lead to noticeable revisions for many nuclear properties given 
by conventional force models. 
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I. INTRODUCTION 



In spite of the great progress attained recently in the construction of the modern reahstic 
iViV-potentials of second generation, based on the concept of improved one- and two-meson 
exchange HTj-Q, a large number of unsolved problems are still left in the field. The majority 
of the problems here is related to description of the short-range part of A^A^-interaction 
and to the quantitative description of few-nucleon systems In particular, one of the 

basic problems is connected with the consistent incorporation of gluon- and quark-exchange 
degrees of freedom and their "matching" with the meson-exchange concept. 

One of the most fundamental difficulties here is how to avoid double counting of the same 
effects in the gluon- and meson-exchange sectors of the unified interaction. Another basic 
problem is the fact the true six-quark microscopic Hamiltonian is presently unknown. While 
some effective three-quark Hamiltonians which include chiral symmetry breaking and con- 



finement [|TO|-0 have been developed to describe the baryon spectra, there is no guarantee 
that the same Hamiltonians can also be applied for six- (and multi) quark systems. 

Recently some fundamental problems have also been found in the consistent de- 
scription of few-nucleon and meson-few-nucleon systems. On the one hand, the numerous 



calculations made in recent years for quark-effects in various few-nucleon observables |T3|jn 
have shown, in general, quite moderate contributions of such effects at low energies and 
momentum transfers [|l^ -[T6[|. On the other hand, however, quite remarkable disagreements 
between the data and the most accurate three- and four-nucleon calculations have been 
found |]5|-^. They probably can be ascribed to an improper treatment of the quark degrees 
of freedom. This follows from the fact that the above mentioned few-nucleon calculations 
do include 3N-force and A-isobar effects together with the most realistic A^A^ interactions, 
i.e. they include, to our current knowledge, all essential contributions. 

Thus, developing a quantitative A^A^-interaction which includes properly both quark- and 
meson exchange effects and is, on the other hand, not so difficult to handle and complicated 
as those A^A^ models derived from multiquark Hamiltonians, is still very topical. Besides, in 
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view of our insufficient knowledge of the accurate six-quark Hamiltonian, the theory must 
be constructed in a manner to avoid features of QCD which presently can not be handled 
with confidence, e.g. in particular the qq interaction, the form of confinement (e.g. linear 
or quadratic) etc. At the same time, however, it is highly desirable to incorporate into the 
interaction model some general reasoning about the preferred symmetry of the six-quark 
system in various NN channels, about the characteristic size of six-quark states etc. In such 
a case the conclusions derived from the model will result mainly from general symmetry 
requirements and general structure of the model etc. rather than some particular choice of 
parameters for the qq interaction, or for the particular law of confinement etc. 

We will demonstrate in the present work that only a few basic assumptions are quite 
sufficient for the derivation of such a hybrid model. On this basis quark effects in the few- 
nucleon physics can be described more reliably. Our consideration is based on the assumption 
(which is common for all hybrid models) that both nucleons merge somehow their quark 
contents at short ranges into different six-quark states dependent on the partial wave, the 
energy and the total spin. While at intermediate and large distances where the nucleons do 
not overlap noticeably with each other the interaction mechanism is governed by one-meson 
exchange (which was just the original Yukawa idea about the origin of strong interaction 



17| , [18|| ). This is a common basis of all hybrid models and one model is distinguished from 



other one by the way of matching of inner quark- and external meson-exchange channels. E.g. 
in hybrid models due to Kisslinger |jl9[ and Simonov , the matching of both channels is 
done at some arbitrarily chosen hypersphere with radius i?, although the matching conditions 
in both models [lT9| , pO| are quite different. 

In a sharp contrast to such hybrid models we prefer to do this matching in the Hilbert 
space of the six-quark states with different symmetries, where every such a state is con- 
structed from single-particle harmonic oscillator (h.o) quark orbits pl|^3|. Thus, in ac- 
cordance to this idea we subdivide a total Hilbert space 7i of the six-quark states on two 
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mutually orthogonal subspaces Tip and Tiq: 

Ti = Hq + Hp. 

In the first subspace Tig we include six-quark states with highest possible spatial sym- 
metry, where one has a maximal overlap of the single-quark orbitals. Whereas the six-quark 
states with a lower spatial symmetry are placed in its orthogonal complementary subspace. 
The most symmetrical states can be shown to have the structure which is rather similar to 
compound states in a spherical or weakly deformed bag. On the other hand the states of 
lower symmetry include a few p-quark orbitals like s^p^, s^p^ etc. and these mixed sym- 
metry states, being projected out onto A^A^-channel (of unexcited nucleons), result into the 



nodal A^A^-radial wavefunctions Accordingly, their structure is analogous to clusterized 
peripheral states. 

On the basis of above considerations and also of other arguments of symmetry character 



we have suggested in the previous works |24-28| a two- component model for baryon-baryon 



interaction with two mutually orthogonal channels. Then, by subsequent exclusion of six- 
quark compound states one comes to an effective one-channel potential model of Moscow 
type, in which a deep A^A^-potential well includes as its eigenstates the most symmetrical 
six-quark states (to be more precise, their projections onto the A^A^-channel). As a result 
of combining two different components into one channel for the effective interaction, the 
orthogonality condition between A^A^ scattering states and localized six-quark states in such 
a model is satisfied automatically due to the hermiticity of the Hamiltonian. In a consistent 
realization of such a program the wavefunctions of the A^A^- relative motion in the "external" , 
i.e. clusterized channel are generally not to be related to the six-quark wavefunctions in the 
"inner" channel. Moreover, it is very likely they should be wavefunctions belonging to quite 
different Hamiltonians. The underlying dynamics of the most symmetric six-quark states 
must be very tightly interrelated to specific chromodynamic effects such as quark- and gluon 
condensates, instantons, breaking chiral invariance etc. whereas the external channel should 
be describable in terms of meson-exchange. 



Thus, these two different channels can hardly be described consistently by a unified 
Hamiltonian (at least, on the up-to-date level of our knowledge for low-energy QCD). The 
price to perform technically this description is a different dynamics in the multiquark- and 
meson-exchange channels^. And hence it is highly desirable to employ for these two compo- 
nents a two-channel model with mutually orthogonal channels. 

This reasoning justifies our model from the physical point of view. Moreover, its real 
success in description of A^A^ partial amplitudes demonstrated in the present paper allows 
to justify the model a posteriori] 

The structure of the work is as follows. In the Section II we present our approach for 
the construction of the two-component hybrid model of A^A^-interaction and its interrela- 
tion to possible dibaryons. In the Section III we present a realization of the generalized 
orthogonality-condition model (GOCM). In Section IV the above GOCM is constructed and 
the structure of the one-channel potential is discussed. Section V is devoted to a quanti- 
tative description of A^A^-phase shifts in the energy range - 350 MeV. We give here also 
the effective-range parameters and a detailed discussion of the structure of the deuteron. In 
Section VI we discuss specific nonconventional interference and tensor mixing between nodal 
and nodeless wavefunctions and the cut-off parameters for the meson-exchange potentials 
at short range. And finally, the main results of the work are summarized in the Section 
VII. In the Appendix we give the formulas for our interaction model in the momentum 
representation. 



"^It should be emphasized that in currently developed models of baryons in which the qq— inter- 



action is described via one- meson exchange [10,21| these one- meson degrees of freedom are nothing 
else but effective degrees of freedom. Thus, these degrees of freedom in multiquark system will be 
somehow different from those in three-quark system. 
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II. A HYBRID MODEL WITH ORTHOGONAL COMPONENTS 



The nucleon-nucleon interaction at large and intermediate distances is well known and 
can be described by meson-exchange potentials . The internal nucleon degrees of freedom 
(quark and gluon ones, if we start from the quark model for nucleon) do not show up in 
this approach. However, when the nucleons come closer than ~1 fm, a transition of the NN 
system into other channels arises, i.e. the internal nucleon degrees of freedom begin to be 
of crucial importance. As a dynamic model, e.g., a six-quark bag model can be used. 

However we still have no full dynamic model describing all possible states of the two- 
nucleon system. Therefore we divide the full Hilbert space Ti (including both nucleonic and 
non-nucleonic degrees of freedom) into two orthogonal subspaces p^-p8|: 

H = Hnn ffi T^eq (1) 

These subspaces must be orthogonal because the dynamics in them is essentially different: 
one of them - Hnn - includes only nucleonic degrees of freedom and meson - exchange 
dynamics, whereas the other, named Ti-Qg, includes 6g-model dynamics (or QCD-inspired 
dynamics). Accordingly to (1) we introduce two mutually orthogonal projection operators 
Pnn and Pgg- It is important that the (unknown) full Hamiltonian of the system does 
not commutate with Pnn and Peg and contains transitions between the A^A^- and the 6q- 
channels. 

If we suppose the existence of a full Hamiltonian H obeying the 6q Schrodinger equation: 

Hip = Eijj 

one can easily obtain, following Feshbach |^9[|, the effective Hamiltonian for iViV-component: 



^pNN = Pnn^ (2) 

PnnHPnn^nn + PNNHPeg[Peg{E - H)PQg]-^ P^gHPNN^NN = E^nn (3a) 

Peqi^NN = (3b) 
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In accordance with these ideas the effective nucleon-nucleon Hamiltonian /ijvtv = 
PjvAf-f^-PjvAf includes only meson-exchange potentials: 

hNr^ = PNr^HPr,N = t + v''^ (4) 

The second term in eq. (3a) is an effective potential which couples A^A^- and 6g-channels 
and will be designated further as vjygN. With these notations, the effective equation with 
orthogonality condition (3b) takes the form: 

{hNN + VNqN)'ipNN = ElpNN (5a) 

Peqi^NN = (5b) 

As an effective wavefunction in A^A^-channel ipNN one can naturally use the resonating 
group ansatz (RGA): 

^NN = A{(P^(P^Xnn)^ 

in which v^at is the nucleon wavefunction and Xnn wavefunction for the A^A^-relative 

motion obeying the orthogonality constraint (3b). 

The approach above formulated can be considered as a general formal framework for the 
hybrid model of the A^A^ interaction. The parameters of vnqn can be determined from the 
underlying six-quark Hamiltonian plip2|, if one assumes that the effective qq interaction 
is the same in 3q- and 6g-systems, or by fitting the A^A^-scattering phase shifts. (A quite 
similar procedure has been used, e.g., in the quark compound bag model due to Simonov 

There are two essential differences in our approach from other hybrid models. These are 
the orthogonality condition in eqs. (5) and the structure of total Hamiltonian. It should 
be emphasized that equations (3) and (5) are fundamentally different. Eqs. (3) is formally 
derived from the full Schrodinger equation by means of the identity transformations. The 
orthogonality condition does not play a role in eq. (3a) except at E=0. On the contrary, 
the eqs. (5) are model equations, which don't involve the full Hamiltonian H. Therefore 
the presence of the orthogonality condition (5b) is absolutely necessary. 
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III. TWO-COMPONENT MODEL IN FRAME OF CONSTITUEUENT QUARK 

MODEL AND MOSCOW POTENTIAL 



To fill the general scheme (5) with a microscopic content it is necessary to use some 
approximation for the full Hamiltonian H . This can be done in the frame of the constituent 



quark model (see e.g. our previous paper [27]). Symmetry considerations allowed to identify 
the subspace Ti^q. It consists out of square integrable functions ip^q describing the lowest 
6g-bag states with maximal spatial symmetry: |s^[6] > for 5- waves and |s^p[51] > for P- 
waves. 

This choice of Tigg can be justified by several independent reasons [p^-|28[|. E.g., recent 



chiral model calculations |^ have shown that the structure of fully symmetric 6g states 



|s^[6] >, in contrast to the mixed symmetry states as |s^p^[42] >, cannot be described by 
the cluster RGM-ansatz and that they are quite similar to the shell-model ground states 
of magic nuclei. However, the most conclusive argument in favor of the s separation of 
6g-states with high and low symmetry arises from our general understanding of quantum 
chromodynamics, where the effective interactions must essentially depend on number and 
type of quarks. Moreover, if we assume, that some effective bosonisation of initial QCD- 
degrees of freedom occurs in the peripheral area of nucleon and thus this bosonic mode is an 
important component of interquark interaction |TD| one can conclude, taking into account 



a highly nonlinear character of such bosonization, chiral meson fields must play a crucial 
role in the dynamics of the six-quark configurations. Such chiral fields should stabilize 
strongly these six-quark components of interaction. Thus, by approximating the Q-space 
Green function [PQq{E — H)Pq^ with one pole a,t E = Egg one gets a separable form for the 
potential vnqn 

VNqN = PNN\H\lpQg > {E - EQq)~^ < 1pQq\H\PNN (6) 

where 

Esq =< ipeqlHl^jjeq > 
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As a projection operator onto the A^A^-channel one can employ a respective operator 
taken from the resonating group method (RGM): 

Pnn = AliJNi^N > < '^n^n\A (7) 

where ipN is the three - quark function of the nucleon, A is the antisymmetrizer, and M is 
the overlap kerneh 

M =< iIjn^n\A\^n^n > (8) 

With this choice of the projection operator Pnn, the model equation (5) becomes a two-body 
effective Schrodinger equation for the orthogonalized relative motion wave function x{R) 

(t^ + V'^^ + 10^^) X = Ex (9a) 

<g\x>=0 (9b) 

in which 

< R|/ /(R) =< ^eq\H\tlJN^pN > (10) 

< R|^ g{R) =< ^Gql^NipN > ■ (11) 
In a good approximation one can take a delta-function for the overlap kernel A/'(R, R') 



21,22 



Ar(R,R')^^(5(R-R0. (12) 

We emphasize once again that eq. (9a) with the orthogonality condition (9b) is not 
equivalent to the full six-quark Schrodinger equation Hip = Eip. Actually we suppose we 
know only individual parts of the full Hamiltonian: 

- subHamiltonian Knn = tR + f acting in the subspace Hnn and describing the 
meson-exchange interaction between unexcitable nucleons, and 
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- other subHamiltonian H^g, describing the lowest states in 6g-bag (in given case H^g = 

Thus the full six-quark Hamiltonian is needed only for determination of coupling between 
the subspaces in eq.(6). In this model the 6g-bag functions ipQg are not eigenfunctions of the 
full Hamiltonian (otherwise [PQg,H] = and v^g^ = 0). Moreover, it is obvious that the 
sum of the projectors P^n (7) and Pgg = J2 |^6g >< ^6g| is not unity in the full six-quark 
space Ti. Therefore, eq. (9) cannot be formally deduced from the full Schrodinger equation 
and the orthogonality condition (9b) proves to be necessary. 

The effective two-nucleon equation (9a) provides the basis for developing the local and 
nonlocal parts of A^A^-interaction models of Moscow type. The main point here is just the 
orthogonality condition (in S- and P- waves), which results in appearance of nodes in A^A^- 
scattering wave functions, the positions of the nodes being do not depend on energy (at least 
up to laboratory energies E^n ~1 GeV). The term vnqn provides an additional attractive 
interaction ai E < E^g. It has been shown in previous papers p^^8|, that the phase shifts 
and nodal behavior of wave functions typical for eq. (9) are well reproduced by a deep local 
attractive potential with an extra bound state and the respective orthogonality condition 
constraint. So, from this point of view, the A^A^-interaction model, known today as Moscow 
potential, is the simplest local model which ensures the orthogonality between the scattering 
wave functions and the most symmetric 6q states \s^[6] > projected onto the A^A^-channel. 
However the situation for P-waves turns out to be already different. Attempts to achieve 
a satisfactory description of the phase shifts by using a local attractive potential failed for 



these partial waves |25|. Therefore, one needs to use the general orthogonality condition 



model (GOCM) presented here. 



IV. STRUCTURE OF THE POTENTIAL 



Here we give the full version of the A^A^ potential model with the additional orthogonality 
condition in S- and P-waves. The potential is an effective one- component approximation to 
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the two-component model, described in the previous section. Actually we have replaced the 
nonlocal term V]\fqN (attractive at low energies) in eq. (9) by an additional local attractive 
well. 

The total interaction is however highly nonlocal due to the presence of the S- and P-wave 
projection operators which are employed in order to take into account the orthogonality con- 
dition (9b). As a result we do not require locality, this means we have a weaker interrelation 
between the orthogonality condition and the form of the attractive well. This decoupling of 
the attractive potential from the orthogonality condition improves essentially the approach. 
In particular, the quality of the fits for P-waves gets more accurate than in the old-fashioned 
Moscow model with eigenprojection [p^-p6[|. Besides the matrix eigenstate projector in cou- 



pled ^Si — ^Di channels, as was demonstrated in our previous paper can be replaced 
quite accurately by a scalar one-channel projector. 

In order to use a potential with the orthogonality conditions in few-body calculations, 
one has to add the projection operator with a very large positive coupling constant to the 
local part of the potential, in all S- and P partial waves P? . 



For the sake of uniformity and convenience we include similar separable terms, but with 
finite coupling constants, also in some other partial waves {D and F). These terms replace 
the standard spin-orbital part of the interaction (for even-parity waves) and reduce partially 
a strong attraction due to the central part of the local potential. In fact, these separable 
terms imitate a short-range repulsion generated by a;-meson exchange^. We include also the 
tensor interaction which couples partial waves with angular momenta / and / ± 2. It can 
be quite accurately described by a truncated OPE-potential in all partial waves with the 
channel coupling being determined by a truncation parameter. 

In the present version of Moscow potential we have replaced the Gaussian form of the 



^It should be emphasized here that the w-exchange terms in traditional meson-exchange models 
are highly nonlocal due to form factors and energy- and momentum dependence. 
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central potential which have been used in all previous versions of the model p^- p6| , |3T| by an 
exponential one. We have found the exponential form gives a more satisfactory description 
of the phase shifts, in particularly for the ^Si — '^Di-channel (see also the refs. ||32||). 
Thus, the model potential consists out of three parts: 

^^iV7V = < + t^''"''' + t^^^P (13) 

where the local exponent well v''^^ depends on the channel spin and parity: 

VMir) = Voexp{-(5r) + {sl)V^' exp{-f3,r). (14) 

In the state-dependent separable part 

v'^'P = X\if X ifl (15) 

a Gaussian form factor < r\(p >= (p{r) is used: 

<^(r)=iVr^+iexp(^-i(^)'^ (16) 

with normalization condition / (p'^dr = 1. The integer / labels the partial waves. 

For the one-pion-exchange part of the potential the standard form with a dipole form 
factor is chosen: 



'°^^W - iaM)^^ (17) 



With such a form factor choice the OPE tensor potential vanishes at the origin as it 
should. In the coordinate representation the OPE-potential has the form: 

f2 



^p)i.™(/,_,(r)(<ria2)+/r(r)Si,) (18) 

where the tensor operator 

A (o-ir)(a2r) (0-10-2) .^o^ 
^12 - (19j 

and 
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fcii^) = (exp(— x) — exp(— q;x))/x — (a^ — l)a/2 exp(— ax); 



(20) 



/r(r) = exp(— x)/x(l + 3/x + 3/x^) — a^exp(— ax)/(ax)(l + 3/ (ax) + 3/ (ax)^) 



(a^ — l)a/2 exp(— ax)(l + l/(ax)); 



(21) 



X = mr\ 



a = A/m. 



(22) 



We use here the averaged pion mass m = (m,ro + '2m.„^)/3 and the averaged value of 
pion-nucleon couphng constant /^/(47r) = 0.075 as we don't wish to deal with the difference 
between np and pp isovector phase shifts in the present work. 

Thus only three free parameters Vq, f3, and a are left for the local part of interaction 
for each combination of spin and parity in addition to two parameters tq and A of the 
separable term in each channel. It should be noted that only some of the values tq and A are 
independent free parameters (for D- and F-waves). Values of A for S- and P- waves must 
go to infinity (in real calculations the value of A ~ 10^ — 10^ MeV is quite enough). Values 
for ro for these channels are related to the local attractive well (for the local potential, the 
requirement of the best approximation for eigen bound state by Gaussian (16) defines ro 
uniquely). Thus, we have totally 32 parameters of the potential (and the value of nNN 
coupling constant) giving a very good description of all — partial waves (except of 
some high / channels) in the wide energy range — 400 MeV. The number of parameters 
almost coincides with that for most recent version of the Nijmegen N — N potential [Q. The 
parameters for the present version of our A^A^ potential are given in Tables I-II. 

V. DESCRIPTION OF PHASE SHIFTS AND DEUTERON STRUCTURE 

The potential parameters as given in Tables I-II were determined by fitting the Nijmegen 



phase shifts (PWA93) §. In Figs. 1-3 the recent SAID phase shifts (SP97) |5§ are also 
presented for comparison. As can be seen from the Figures, some discrepancy between the 
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results of both phase shift analyses (PSA) exist, especially for some partial phase shifts. 
With applications to few-nucleon problems in mind we tried to reproduce with maximal 
accuracy the ^Sq and ^Si — ^Di phase shifts and the values of the scattering length and the 
effective range. 



A. Singlet partial wave channels 

The description of singlet n — p phase shifts for both even- and odd parities is illustrated 
on Fig.l. It is evident from the Figure that the quality of fit to the data of recent phase 
shift analysis is quite good, especially for the Nijmegen PSA-results. E.g. the fits in ^5*0 and 
^Pi channels are almost perfect. The quality of fits can be estimated quantitatively from 
the Table IV for these channels. The average deviations for all singlet channels are only 
0.1 — 0.2% excepting the ^G'4-channel where the discrepancy with PSA-data is largest and 
around 1%. 

Also there is a problem in precise description of the singlet effective range tq (see the 
Table III). We used as "experimental" the value tq presented in compilation of Dumbrajs 
et all of 1983 (see footnote to the Table III). However, in view of the very good agreement 
of our phase shifts with the Nijmegen PSA for ^S'o-channel one could conclude that the 
disagreement for tq should be really much reduced. 



B. The even-parity waves 



The S'-wave potential turns out, as is in the previous versions of Moscow potential ||27|| , 



to be strongly attractive. The Gaussian (15) with the range parameter ro, included in the 
orthogonality condition (9b), is close to the eigenfunction of the ground "forbidden" state 
in the potential. In other words, we obtain for S-waves actually almost a local poten- 
tial. However, for the ^Si — ^Di channel (and also for all triplet coupled channels) we use, 
strictly speaking, non-eigenstate one-channel projector, as in in order to avoid a more 
complicated two-channel eigenprojector. 
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We do not introduce here a spin-orbital potential for even partial waves in an explicit 
form because it cannot be determined by PSA-data for ^5"! — ^Di channel, and the role of 
spin-orbital potential for higher even-parity partial waves is played by the term v^^p. 

It should be kept in mind here that the complete two-channel version of our model 
includes in the proper A'"A'"-channel one-meson exchange interaction terms (in a subspace 
orthogonal to symmetric six-quark compound states). Thus, in the two-channel model, the 
spin-orbit terms should be described by a conventional meson-exchange model. However in 
the effective one-channel model presented here the separable state-dependent spin-orbit in- 
teraction in even-parity channels is inavoidable to compensate partially the strong attractive 
potential in the S-wave. 

The effective range parameters for singlet and triplet S-wecve channels are given in 
Table III. Among all the calculated phase shifts the maximal disagreement with PWA93 
(though not large) is observed for the tensor mixing parameter ei. 

C. The triplet odd-parity waves 

In the accordance to the microscopic quark picture, the orthogonality to the bag-like 
functions |s^p[51] > must be included for all P-waves. However, if we shall look at the 
behaviour of "experimental" P-wave phase shifts at energies up to 500 MeV we shall not 
find any repulsion in the channel, because the corresponding phase shifts are purely 
positive until the energies ~1 GeV. There is no repulsive core for this channel also in the 
majority of the conventional realistic A'"A^ potentials. But if we look to the phase shifts at 
high energies (see Fig. 4) one can observe a repulsion appearing in all three triplet P-waves, 
while ^P2-phase shifts becoming negative only at energies higher 1 GeV. Prom the point of 
view of the constraints imposed by the orthogonality condition, this means the function to 
which the ^P2 scattering function is orthogonal is much more narrow than that for other 
P-wave channels, ^Pq and ^Pi. 

It is interesting that fitting the ^P2-wave at energies up to 350 MeV enables us already 
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to determine the range parameter ro of the projector (see Table II). The inclusion of the 
projector improves appreciably the description of the phase shifts up to 350 MeV. An attempt 
to reproduce the ^P2-phase shift using a purely attractive potential with "extra" bound state 
results in a very deep (~15 GeV) potential and an unsatisfactory quality of the description. 
Besides, such a deep potential is not suitable for description of the ^Pq and ^Pi phase shifts. 
That is why we have waived in the present version from the concept of a local Moscow model 
for P- waves. 

So, for odd partial waves we have a rather small attractive well (~220 MeV) and orthog- 
onality to the non-eigen bound states for local part of potential. It might mean the size of 
six-quark bag in P-waves should be smaller than in ^Si and ^Sq waves. One notices here 
that the range parameters of the projectors for the ^Pq and ^Pi-channels (ro ~ 0.32 fm) are 
almost coinciding with each other. As is seen from Fig. 4, attraction for some odd higher 
partial waves with L — J {^F^ and ^H^) is noticeably deficient in the given version of the 
model. 

Unlike the even partial waves, we used a more conventional local spin-orbit potential for 
odd partial waves (see eq.( 14)) because the usage of the separable spin-orbital form is not 
convenient to describe the splitting of P-phase shifts. 

It would be rather instructive to estimate the averaged relative difference of phase shifts 
predicted by the Moscow model and the recent phase shift analysis [1] using the criterion of 
relative difference or the respective absolute difference measured in radians: 



1=1 



1 ^ 



^PSA 



2 

(23) 



XjSl — ^ X] 1^^31,1 ~ <^J5M ' (^^) 



i=l 



where Sj% j^ and Sjg^ are partial phase shifts in the channels JSl at the energy Ei for the 
Moscow model and Nijmegen phase shift analysis respectively. The Table IV presents the 
values of Erei and Xjsi for all considered JS'Z-channels. It is evident from the Table the 
average deviation of phase shifts predicted by the Moscow model and recent PSA is very 
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small and around 0.2 — 0.4%. It means that the description of A^A^-observables with the 
presented force model should be very good. 



D. Deuteron structure 

The accurate description of the deuteron structure offers an additional strong test for any 
nuclear force model. Many deuteron properties, even in the static limit, depend sensitively 
on the behaviour of the A^A^ force at intermediate and short ranges |T^, especially on the 



Z^-wave contribution. For example, with the first version of the present force model |^ 
we found an impressive agreement with experimental data for all crucial D-wave deuteron 
observables like Qd, As, Ad/ As- But this early model included a node not only in the S-wave 
but also in the D-wave. This extra node in D-wave was a consequence of a very short-range 



truncation of the OPE tensor force which contradicts somehow the microscopic picture 
of the underlying interactions (e.g. according to the wide-spread opinion |^ the OPE tensor 
force cannot penetrate deeply inside the two-nucleon overlap region). 



Hence, in subsequent versions of the model [p^ -|28[], a more soft cut-off factor has been 
employed which resulted in the disappearance of the D-wave node. As an immediate con- 
sequence of the softer truncation in the OPE tensor force the D-wave deuteron observables 
have become close to the values predicted by conventional force models, i.e. the values of rj 
and Qd are a little bit underestimated (see Table V). Nevertheless the node in the 5- wave 
and the strong attractive S-wave potential, tightly related to this, results in a very specific 
interference between S- and D-wave components and a specific character of tensor mixing 
(see Section VI). 

The values of deuteron observables for three versions of our force model are presented 
in Table V while the pattern of the deuteron wave functions is displayed on Fig. 5. One 
can see on the Figure the short-range maximum in the D-wave almost disappears for the 
current version of the force model while this maximum in the S-wave gets rather reduced. 
It is interesting to note the D-wave amplitude in the current version of the model (solid 
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line) is a little bit lower than in the previous versions (dashed and dot-dashed lines) due to 
a smaller value of the derivative of the D-wave component near the S-wave node (~ 0.53 
fm). While the asymptotic behaviour of the 5*- wave looks almost perfect (see values of As 
in Table V). 

Thus we can conclude from the deuteron results presented in this Section that the short- 
range part of the tensor force needs to be a bit improved. Careful inspection of the Table V 
shows unambiguously the general good agreement for the deuteron parameters found with 
the sharply different force models such as Nijmegen and Moscow potentials. The values for 
the deuteron observables are a result of some general properties (like OPE tail) and of the 
A^A^ phase-shifts used for fitting only and essentially do not depend on the details of the 
force at short ranges^. 

VI. SPECIFIC INTERFERENCE BETWEEN TENSOR AND CENTRAL 
FORCES AND THE ttNN COUPLING CONSTANT IN THE MOSCOW FORCE 

MODEL 

We included this specialized Section to the present work in order to emphasize a specific 
character of interference between tensor and central forces in the Moscow force model. This 
interference will be shown below to be very advantageous in some aspects as compared 
to the traditional force models. The main difference between our and traditional models 
as concerned to wave function form is the nodal character of the S-wave deuteron and 
scattering wavefunctions and the practically nodeless character of the D-wave functions.^ 



Certainly this conclusion may be invalid for non-static, e.g. energy-dependent or multicomponent 
force models. 

^The very small inner maximum in the D-state wavefunction in the present version (see the 
solid lines on the Figs. 6a-6b) can be ignored in any calculation if we do not consider very high 
momentum transfer. 
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We will show here that the specific tensor mixing between the 5'-wave state with a node and 
the almost nodeless D-wave state results in a remarkably different ei-behavior. We compare 
the D-wave observables with the results of traditional models. 

First of all we emphasize here that the best fit for A^A^ phase shifts is attained in our 
case with a very reasonable value for the OPE cut-off parameter K^ip = 0.78 GeV (we used 
here the dipole form factor), see eq.(17). This soft cut-off parameter is in nice agreement 
with both experimental results and with all theoretical estimations made in vr-N dynamics 
|33|^7|. It should be contrasted with a statement formulated in |3^, p. 232] for traditional 
OBE-force model: "...a value of 1.3 GeV is a lower limit for A^r". The conventional OBEP 
model with A = 0.78 GeV gives the extremely low values for = 0.238 fm^, the ratio 



D/S = 0.0233 and Pd = 2.4% which should be compared to the respective values for 
our force model (see Table V). 

In despite of the "soft" value of A, the D-wave deuteron properties in our model (see 
Table V in Sect. V) are in a rather good agreement with the experimental data, being 
remarkably better than the respective predictions of the traditional force models with the 
same A value. We note, in passing, that the harder truncation with A ~ 1.3 ^ 1.7 GeV 
is usually taken in the traditional force model just in order to fit reasonably the deuteron 
properties and the tensor mixing parameter (see below). 

The second important point in the story is related to the mixing parameter ^i. In fact, 
in order to reach a reasonable agreement with the recent phase shift analysis data for the 
^i-mixing parameter |]I|,|33| the A- value must be taken also around 1.5 1.7 GeV ||3^ (see 
Fig. 6) while the same agreement with the experimental ei is reached in our model using a 
much more soft A = 0.78 GeV. This sharp difference from the traditional force models can 
be ascribed to a different character of mixing between 5*- and Z^-waves in our model. 

Some additional confirmation comes from the value of vrA^A^-coupling constant obtained 
in our model. We choose here the Nijmegen force model as a good representative of the 
traditional A^A^ potentials (see, the Table V). Two above models include practically the 
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same values for the ttA^A^ (charged) couphng constants {f^^N = 0.075 in our case^ and 
fw±NN ~ 0.0748 for Nijmegen potential). The latter fact is very important because the D- 
wave characteristics are directly related to the ttNN coupling constant. In this respect our 
model appears to corroborate the smaller value of qI^n — 13.60 advocated by the Nijmegen 
group 1T|§. 

The two nice features of our model discussed above, i.e. the soft cut-off parameter A and 
low value of irNN coupling constant, are in agreement with modern trends and lend strong 
support to our model. 



VII. CONCLUSION 

The force model presented in this paper differs in a few important aspects from tra- 
ditional A^A^ interaction models currently in use. First of all the Moscow two-component 
model includes two mutually orthogonal quark- and meson-exchange channels. This channel 
orthogonality leads to many differences from the traditional force models. In particular it 
requires a node in low partial waves with the node position almost independent on the rela- 
tive energy in a wide energy range (< 1 GeV). The nodal behaviour of wave functions is also 
preserved for the one-channel model presented here. The node in the A^A^ wave functions 
results in an enhancement of high momentum components and a strong increase of the aver- 
age kinetic energy in the deuteron and in all few-nucleon systems. This increase of the inner 
kinetic energy leads to significant enhancement of higher angular momentum components 
in nuclei and nuclear matter and also for many particular nuclear processes like vr- 



meson absorption and scattering in the A-resonance region etc. This strong enhancement of 
high-momentum components in A^ — A^ system as compared to any traditional N — N force 



model may be seen e.g. in hard bremsstrahlung process pp pp'j [39| at Ep = 300 MeV 



°The value corresponds to the charged coupling constant because we considered first of all the pn 
scattering phase shifts. 
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and higher at small forward and backward angles 9^ of 7-emission. To make the compar- 



ison with traditional repulsive core models most unambiguously the authors |^ did their 
bremsstrahlung calculations with both the Moscow model (in its previous version [^) and 
its exact phase-shift equivalent supersymmetrical partner. Thus, such a comparison removes 
any questions on the possible on-shell origin of disagreements observed. 

Redistribution of higher partial waves along Jacoby coordinates leads e.g. to a notice- 
able enhancement of the P-wave attraction for N + d and + 2a systems [^,0. The 
long-standing puzzle of the analyzing power Ay in low energy N + d scattering is explained 
by insufficient attraction in just the N — d relative motion P-wave The apparent 

discrepancies for n + elastic and n + '^He d + d rearrangement low-energy scattering 
observed recently ||^ appear to have to be explained also by insufficient attraction in the 
n + '^H(^He) P-wave Such enhancement of higher partial wave contributions to near- 

threshold- and low-energy processes in few-nucleon and few-cluster physics when replacing 
the deep Moscow-type potential (including extra bound states) with its SUSY-partner po- 
tential - which is exact phase-shift equivalent - is a sequence of some very general algebraic 
properties of kinetic energy operator in different coordinate systems and is disconnected at 
all to any small variations in the on-shell properties of various N — N potential models of 
current use. 

The second crucial point in the development of Moscow A^A^ force model is the important 
role of the six-quark components with maximal possible symmetry. We showed recently that 
the coupling of the meson-exchange NN channel to the six-quark component can be strong 
enough to represent adequately the intermediate-range A^A^ attraction. In turn, this fact 
leads to quite remarkable contributions of such six-quark configurations in nuclear bound 
and low-excited states. If so, it may require some strong revision for many nuclear properties 
as given by traditional force models (e.g. the meson-exchange current contributions). Thus 
the strongest test for the new model may offer few-nucleon calculations for the analyzing 
power Ay in the n + d and p + d low-energy scattering, for the analyzing power Ay in p + d 
radiative capture reaction and for the p + d intermediate energy elastic scattering cross 
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sections (so called the Sagara puzzle Hence the careful comparison of the predictions 

for few-nucleon systems using the Moscow force and more traditional NN interactions may 
be extremely interesting. 
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APPENDIX: MOSCOW POTENTIAL IN MOMENTUM SPACE 



The X-matrix defined as 

2iMqk = (Al) 

( M is the reduced mass while g is a hnear momentum) obeys the partial-wave Lippmann- 
Schwinger equation: 

where P means the principal value integral. The elements of matrix V in eq.(2) are equal 
to the partial- wave momentum-space potential in IS J basis (up to factor 1/A.ti): 

Vn(q',q) = ^< (l'S)J\V(q-c^)\(lS)J>, (A3) 

where V{q) is related to V{r) by a standard Fourier transformation: 

T/(q) = y e-^(^'^V(r)dr (A4) 

Here we give explicit formulas for all terms of the present version of Moscow potential 
Viii{q',q) in momentum space (in MeV~^). 



1. Local part of Moscow potential VJ^! 



loc 



(A5) 



where the parameters (3 and /3i are given in MeV: 

$ — P ■ he, (3i — Pi ■ he 

Fi is being the derivative of the second kind Legendre function: 

d ^, . „ , , 1 dzPi{z) 
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2. Separable terms of the potential 

In momentum space the separable terms with Gaussian form factors (15,16) have the 
same form as in the coordinate space: 

vrM)-h'>^l<pM)<pM'l (A7) 

where 

Here the normahzation condition / (pf{q)q'^dq = 1 is assumed and the factor 7r/2 is related 
to the integration measure used in eq.(A2), and fo is given in MeV"-*^: 

^0 = ro/{hc). 

3. The OPE potential with dipole truncation 

For the sake of reader's convenience we give also the known formulas for OPE matrix 
elements. 

a) The central part of OPE potential: 
Here and below 

q' + q'' + ml q' + q" + A' 

2qq' ' 2qq' ' 

b) The tensor part of OPE potential for triplet uncoupled channels with / = J: 

V ten )JJ\H,H) 3 47rm2\ qq' 2 J + 1 2 J + 1 -^+7 ^ ' 

where the function Gi is introduced as follows: 

Gi{q. q') = Qi{x) - Qi{y) - {y - x)Fi{y) (A12) 
and X and y are defined by eq.(AlO) 
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c) The tensor part of OPE potential for coupled channels with I — J ±1: 

(A13) 



(A14) 



3 //(j + i) r 

(Ken W+i(5,5) = ^-4^;;^^7qr3- (A15) 



{V,^r)j+.,j-i{ci. q') = (K?/''^).-i,.+i(g', q) (A16) 
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TABLES 



TABLE I. Parameters of local part of the potential 



spin 


singlet 


singlet 


triplet 


triplet 


parity 


even 


odd 


even 


odd 


a 


6.08671 


6.08671 


6.08671 


4.3160 




-4346.19 


-1767.26 


-4567.12 


-223.63 




3.49366 


2.84152 


3.81272 


2.4959 










-591.1 


Pi 








3.4688 




TABLE IL Parameters of projectors and 


separable parts of the potential 




State 




A, MeV 




ro, fm 


'So 




00 




0.3943 


'Pi 




00 




0.5550 


'D2 




107.2 




0.4527 


'Fs 




182.6 




0.5191 


'Si 




00 




0.3737 






161.2 




0.4695 






588.2 




0.3572 


'G, 




2.74 




0.8077 


'Po 




00 




0.3209 


'Pi 




00 




0.3226 


3p 




00 




0.1632 






5.447 




0.6221 
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TABLE III. Effective-range parameters for the potential variant given in Tables 1-2 



a, fm ro, fm 





theory 


experiment 




theory 




experiment 


triplet 


5.422 


5.419(7)^ 




1.754 




1.754(8)^^ 


singlet 


23.74 


-23.748(10)'^ 




2.66 




2.75(5)^ 


^S. Klarsfeld, J. Martorell, and D.W.I. Sprung, J.Phys. 


G: Nucl.Phys. 10, 165 (1984) 




''O. Dumbrajs et al, Nucl.Phys. B216, 277 (1983) 














TABLE IV. Accuracy of fitting of phase shifts 






channel 


'So 


'Pi 


'D2 


'Fs 


'G, 


'H, 


's. 






0.007504 


0.000527 


0.002113 


0.000197 


0.012030 


0.001802 


0.004524 


0.000867 


per point^ 


0.005282 


0.001731 


0.000627 


0.000008 


0.000200 


0.000055 


0.005595 


0.003922 


channel 


£i 






'Gs 


£3 


'G, 






^rel 


0.006816 


0.000022 


0.034310 


0.027256 


0.023088 


0.001451 


0.000184 


0.000455 


per point 


0.000843 


0.000007 


0.000856 


0.004545 


0.013887 


0.000029 


0.000186 


0.003596 


channel 


3p 










'H, 


£4 




^rel 


0.003874 


0.010679 


0.007830 


0.017762 


0.007590 


0.021424 


0.010617 


0.026970 


per point 


0.018245 


0.000123 


0.000693 


0.001413 


0.000279 


0.000004 


0.000042 


0.000155 


spot XPSA „ 
1 V^W^ ("jSl.k "jSl,k\^ 

°JSl,k 

- S^i^f (in radians) 
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TABLE V. Deuteron parameters for conventional and Moscow NN potentials 



model 


Ed (MeV) 


Pd (%) 


rm (fm) 


Qd (fm2) 




As (fm-V2) 


D/S 


^loop 


RSC^ 


2.22461 


6.47 


1.957 


0.2796 


0.8429 


0.8773 


0.0262 




1 N 1 J 111 C/ KJ 




5 754 


1.966 


0.2706 


0.8429 


n 8844 


n 09594 




Moscow 86^ 


2.22444 


6.57 


1.966 


0.2862 


0.8422 


0.8838 


0.0268 


0.53 


Moscow 98^^^ 


2.22440 


5.75 


1.954 


0.2708 


0.8470 


0.8746 


0.0259 


0.30 


presenf^ 


2.22456 


5.65 


1.967 


0.2731 


0.8476 


0.8845 


0.0255 


0.08 


experiment 


2.224575(9) 




1.9660(68) 0.2859(3) 0.857406(1) 


0.8846(16) 


0.0256(1)^ 





'^Dioop is the relative amplitude of the D-wave maxima, i.e the absolute value of the ratio of the 

first and second maximum of the deuteron D-component. 

^The value h'^/2m = 41.47 MeV-fm^ has been used (m = 938.978 MeV). 

^The value h'^/2m = 41.47107 MeV-fm^ is used (m = 938.918 MeV). 

'^Unfortunately, in our previous work [27| only rounded values for potential parameters are given 
in the Table III. The deuteron parameters cited in |2^ (for variant B) do not correspond to the 
rounded potential parameters cited in Table III of ref. [E^. We thank Dr. S.B. Dubovichenko who 



has attracted our attention to this disagreement and give here the exact values for variant (B) of 
Ref. [27]: Vq = -1329.18 MeV, r/ = 2.2959 fm~^, a = 1.8835 fm^^ 
''The present value is taken from |41| 
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Figure captions 

Fig.l. The comparison of spin-singlet phase shifts for the present version of Moscow 
A^A^ potential with the data of the recent energy-dependent phase-shift analyses: PWA93 
g (circles) and SA1D97 (triangles). 

Fig. 2. The spin-triplet even-parity phase shifts for the present version of Moscow A^A^ 
potential. The data of the energy-dependent phase-shift analyses are: PWA93 |l[| (circles) 
and SAID97 H (triangles). 

Fig. 3. The spin-triplet odd-parity phase shifts for the present version of Moscow A^A^ 
potential. The data of the energy-dependent phase-shift analyses are: PWA93 |l|] (circles) 
and SAID97 Q (triangles). 

Fig. 3. (Continued.) 

Fig. 4. The spin-triplet P-wave phase shifts in a wider energy region: the data of energy- 
dependent phase-shift analysis SAID97 |^ (solid lines) and predictions for present version 
of Moscow A^A^ potential (dashed lines). 

Fig. 5a. The deuteron S-wave and D-wave functions for present (solid lines) and previous 



versions ( [^] - dashed lines and [p7| , variant B] - dot-dashed lines) of the A^A^ Moscow-type 
potential. The deuteron wave functions calculated with the RSC potential (dotted lines) are 
shown for comparison. 

Fig. 5b Short- distance zoom of Fig. 5a (see caption to the Fig. 5a). 

Fig. 6. The energy dependence of the mixing parameter ei for different values of cut-off 
parameter A corresponding to conventional (dashed lines) and present (solid lines) force 
models. 
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